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CASE  STUDIES  OF  FATIGUE  CRACK 
GROWTH  USING  AN  IMPROVED  MICRO-LIGAMENT 
INSTABILITY  MODEL 


1.  INTRODUCTION 

A  current  emphasis  of  the  crack  propagation  phase  of  fatigue  failure  is  justified  in  several 
ways.  U  is  seen  to  underlie  estimates  of  air-frame  structural  integrity  and  damage  tolerance. 
The  fatigue  stress  allowables  of  welded  structures,  such  as  offshore  oil  recovery  platforms  and 
pipelines,  appear  to  be  largely  governed  by  crack  propagation  life.  The  ASME  pressure  vessel 
code  utilizes  crack  propagation  rate  data  for  assessing  the  hazard  of  positive  NDT  indications. 
Even  in  aircraft  jet  engine  parts  retired  on  reaching  a  given  probability  of  small-crack  initiation, 
the  remaining  propagation  life  is  of  interest  in  assessing  the  hazard  of  emergency  excess  hours 
without  inspection,  or  the  risk  of  continued  use  of  overaged  parts  which  pass  inspection. 

Fatigue  propagation  life  of  a  structure  is  best  measured  by  direct  simulation  of  its  loading 
spectrum.  However,  it  is  less  costly  to  estimate  this  from  steady-state  data,  modified  by  suit¬ 
able  overload  retardation  models.  Even  here,  the  cost  and  time  required  to  develop  sufficient 
data  can  be  burdensome.  A  goal  of  propagation  modeling  is  to  reduce  this  burden.  With  physi¬ 
cally  based  models,  an  additional  benefit  could  be  an  enhanced  understanding  of  microstructural 
and  environmental  effects. 

As  with  previous  NRL  reports  on  this  endeavor  [1,2,3, 4),  this  one  puts  forth  modeling 
refinements  that  evolved  in  attempting  additional  data  analyses.  In  the  present  set,  29  combi¬ 
nations  of  composition  and  heat  treatment  of  alloys,  comprising  116  variations  of  conditions  of 

load  ratio,  frequency,  and  environment,  are  analyzed.  Some  75%  of  these  materials  are  newly 
Manuscript  submitted  December  10,  1979 
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reported  here,  the  rest  re-analyzed  to  provide  a  more  complete  set.  Except  for  one  m^jor 
modification  affecting  stage  II  R-effects,  the  present  model  is  essentially  as  that  detailed  in  Ref. 
[4].  Accordingly  this  report  presents  only  essential'  steps  of  the  derivation,  which  has  been 
reduced  to  a  simpler,  more  direct  form  than  before. 

It  should  be  stated  at  the  outset  that  the  development  to  follow  here  is  not  a  rigorous 
analysis  of  the  elastic-plastic  crack  tip  field  nor  of  the  micro-separation  processes  which  reside 
in  it.  It  is  rather  a  refined  dimensional  analysis,  combining  simple  limiting  solutions  of  these 
problems  in  a  manner  that  permits  correlation  of  a  large  body  of  experimental  data.  Success  at 
correlation  is  only  circumstantial  evidence  of  proper  combinations.  However,  while  awaiting 
more  rigorous  solutions,  this  approach  is  thought  to  provide  useful  insights,  and  even  predictive 
capabilities,  otherwise  unavailable. 

2.  MODELING  EQUATIONS 

This  model  explores  the  proposition  that  stability  of  crack  propagation  requires,  more 
essentially  than  other  things,  maintenance  of  a  state  of  constant-load  stability  within  a  material- 
specific  region  very  close  to  the  crack  tip.  For  simplicity  of  modeling,  it  considers  the  material 
in  this  region  to  be  formed  into  tiny  cylindrical  tensile  ligaments,  which  are  extended  by  the 
external  loading  as  well  as  by  the  consequent  stable  propagation  of  the  crack.  As  any  specimen 
subjected  to  tensile  strain,  strain-hardening  enhances  stability,  while  Poisson  contraction 
reduces  it.  Two  other  agents  of  instability  are  germane.  Firstly,  when  active  straining  ceases  at 
the  crack  tip  during  the  dwell  time  of  cyclic  loading,  stress  relaxation  can  reduce  ligament  load 
bearing  capacity,  hence  stability.  This  is  a  volume  effect  which,  except  for  bulk  diffusion 
processes,  should  be  independent  of  the  environment.  Secondly,  because  the  ligament  is  so 
small,  environmental  attack  can  cause  a  significant  reduction  of  its  effective  area,  hence  load 
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capacity.  Such  attack  can  be  considered,  in  the  manner  of  corrosion  processes,  as  a  surface 
effect  limited  by  exposure  time  up  to  a  saturation  condition.  It  may  be  limited,  alternatively,  by 
insufficient  surface-disturbing  plastic  strain  due  to  crack  loading  and  propagation. 

2.1  Limited  Tensile  Instability 

Consider  a  simple  tensile  specimen  of  sectional  area  A,  subjected  to  load  P,  for  which  true 
stress  <r  is  defined 

a  —  PI  A ;  P  —  a  A ;  dP  “  adA  +  Ada  (1) 

A  condition  verging  on  tensile  instability  is  dP  *  0,  at  a  load  maximum,  for  which  Eq.  (1) 

becomes 


adA  —  —  Ada  (2) 

To  replace  the  total  differentials  of  area  and  stress,  on  a  ligament  of  radius  rT  -  dTl 2,  the  fol¬ 
lowing  partials  are  germane: 


dA„  •*  —  2v  it  rfdc,  dAs  **  —  2tr  rTdrT  (3) 

do>  ■»  Ode,  d<rm  “  —  ma  dt/t  (4) 

Here  v  is  the  Poisson  ratio,  t  is  strain,  9  is  the  strain  hardening  rate  (da/de)%  m  is  the  stress 

relaxation  exponent  and  t  is  time.  Combining  these  in  Eq.  (2)  and  solving  for  de 


d  €  — 


(5) 


2.2  Crack-Tip  Straining 

In  Eq.  (5),  de  may  be  regarded  as  the  strain  differential  required  to  maintain  the  constant 
load  condition,  dP  —  0.  An  expression  for  how  such  strain  derives  from  the  loading  and 
growth  of  the  crack  is  needed.  This  model  employs  the  simplest  of  crack  tip  plasticity  models, 
the  analogue  of  the  linear  elastic  stress  field,  providing  an  inverse  half  power  singularity 
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strength  near  the  crack  tip.  It  is  employed  provisionally,  as,  for  materials  of  ordinary  strain  har¬ 
dening  exponents,  several  analyses  have  indicated  stronger  than  half  power  singularities.  How¬ 
ever,  use  of  this  weakest  possible  singularity  seems  to  provide  closer  data  correlations,  in  con¬ 
junction  with  other  assumptions  of  this  model  algorithm.  The  analogue  used  is 

<ry  =  K,  /  yJ2nr  ;  e  —  K/ZJlnr  E  (6) 

whence  differentiating 

dt  -  [r-'/2  dK-j-  r~i!2  dr  /JR  E  (7) 

where  K,  is  the  opening  mode  stress  intensity  factor,  r  is  distance  ahead  of  the  crack,  and  E  is 
Young’s  modulus.  Let  r  —  dT,  dr  —  —  da  where  a  is  crack  length,  then 


dt  —  dtL  +  — ^ —  da  (8) 

4rr 

Here  dtL  designates  a  strain  differential  due  to  crack  loading,  while  eG  a  strain  associated  with 
the  strain  gradient  dt/dr ,  hence  —  dt/da ,  at  the  point  r  —  dT.  Since  we  need  be  concerned  with 
maintaining  stability  only  after  the  loading  has  ceased  in  the  fatigue  cycle,  we  neglect  dtL  in 
substituting  Eq.  (8)  in  Eq.  (5),  yielding 


|8  drT  -f  4 m  ry 


dt  1  (  6  .  1  ' 

7 1  ec|7_I 


2.3  Specifics  for  the  Growth  Rate  Equation 


The  general  crack  growth  relationship,  Eq.  (9),  shows  that  the  principal  impediment  to 
crack  growth  is  the  strain  hardening  rate  of  the  material.  Ordinary  and  cylic  strain  hardening 
characteristics  differ  markedly.  Both  are  involved  in  some  way  at  the  crack  tip  in  fatigue. 
Cyclic  loading  produces  cyclic  straining  at  the  crack  tip.  Yet  if  there  is  substantial  growth  dur¬ 
ing  the  period  of  sustained  load,  the  crack  may  invade  material  of  virgin  characteristics, 
unaffected  by  the  cycling.  After  trials  of  many  correlation  algorithms,  the  option  chosen  is  to 
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divide  the  effects  decisively  on  the  basis  of  limiting  cases  of  crack  propagation,  namely:  to  asso¬ 
ciate  cyclic  strain  hardening,  a  function  of  the  strain  excursion,  with  the  impedance  of  stress- 
relaxation  induced  growth;  and  to  associate  monotonic  strain  hardening,  a  function  of  the  max¬ 
imum  strain,  with  impedance  of  environmentally  induced  growth.  This  defines  two  distinct 
forms  of  the  square-bracketed  last  term  of  Eq.  (9),  called  a  growth  rate  factor  G 


G,  = 


and 


G,  = 


ft £  >/3 

€1 

o-,  2 

02p  V3 

€  2 

<t2  2 

-i 


(10) 


(11) 


In  the  foregoing,  9]r  and  02r  are  plastic-only  strain  hardening  rates,  =*=  0~*  -  E~\ 
corresponding  to  the  plastic  strain  dependency  of  the  stress  relaxation  and  the  work  hardening 
processes.  Here  the  monotonic  stress  <T|  is  as  measured  from  a  zero  stress  origin,  and  con¬ 
verted  to  true  stress.  However  the  cyclic  stress  is  the  positive  excursion  of  a  cyclic  loop  in  bal¬ 
ance  with  its  negative  excursion,  as  disposed  in  a  centroid-neutralized  equilibrium  position,  as 
jusitifed  in  Ref.  [51.  The  term  \/3/2,  slightly  less  than  unity,  can  be  justified  as  a  plane  strain 
triaxiality  effect  on  the  instability  condition  [6].  It  has  the  effect  of  increasing  the  point  of  ter¬ 
minal  instability  in  a  way  helpful  in  correlating  materials  with  very  flat  monotonic  stress-strain 
curves.  With  this  growth  partition  hypothesis,  Eq.  (9)  can  be  expressed  as 


da  -  4  rT  m  G2  dt/t  +  8G,  drT  (12A) 

It  is  helpful  at  this  point  to  express  the  crack  growth  as  a  finite  increment  in  a  cycle  of 
crack  loading,  and  corresponding  ligament  straining.  Consider  loading  from  A'mln  to  A'nm  in  an 
interval  tL%  and  thence  holding  on  Kniax  for  an  interval  tH.  With  this,  Eq.  (12)  can  be 
expressed,  assuming  the  change  in  rT  to  be  relatively  small 


A<z  —  A  rT  m  ln(l  +  t2/f tj) G 2  +  8  A  /y  G | 


(12B) 
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or  abbreviated 


Ao  -  S\  G:  +  / 1  G,  (120 

Here  incremental  steps  approximate  the  differentials.  A  summed  or  averaged  growth  rate  factor 

may  be  defined  as 


GRF  -  A«//,  -  G2  +  ^  (7, 
J2 


(13) 


In  this  paper  we  have  used  a  single  value  for  the  loading  wave  form  factor,  .'n(l  +  f/////) 
—  0.3  (vs  0.4  in  Ref.  [4]).  This  value  is  used  to  fit  data  for  both  sinsoidal  and  triangular  wave 
forms  used  in  the  FCGR  tests.  With  this  value,  Eq.  (13)  can  be  restated 


GRF  = 


A  a 


1 .2  mrT 


G>  + 


A  T  j 


C?  |/0. 1 5  m 


(14) 


2.4  Curves  of  Time-  and  of  Strain-Limited  Surface  Attack 

The  two  ways  of  limiting  surface  attack  in  the  fatigue  cycle  were  noted  earlier.  The  effect 
of  varying  these  limits  is  displayed  by  two  families  of  curves.  One  is  a  family  of  curves  for  a 
parametric  set  of  constant  values  of  />-normalized  surface  attack  intrusion,  (Arr/rr);  the  second 
is  a  set  for  such  intrusion  relative  to  the  total  plastic  tensile  strain  on  the  rfrligament  during 
one  loading  cycle,  (Arr/r,/A«J.  Identifying  the  family  members  in  best  fit  to  the  CFCG  data, 
sets  the  value  of  the  fitting  parameters,  other  than  rr. 


For  the  first  family,  Eq.  (14)  is  restated 

GRF  ”  G2  +  (Arr//-rl  CJ i/0. 1 5  m 

or  for  the  parametric  N-family 

G(N)  -  G2+  (21  G|/0.15  m 


The  crack  growth  rate  is  thus 


1.2  m  rT  G(N) 


(15) 

(16) 

(17) 
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For  the  second  family,  Eq.  (14)  is  restated 

GRF  -  G2  +  Ac  —S~—  Gi /0.15  m  (18) 

Ac 

Here  the  Ac  term  is  intended  to  include  plastic  strain  due  to  both  the  loading  and  propagation 
of  the  crack,  as  in  Eq.  (8),  approximated  here  in  incremental  form,  with  Eq.  (14)  substituted 
for  A  a 

Ac  —  Ae^  4-  0.3  m  e0  GRF  (19) 

The  loading  (plastic)  strain  \eL  is  associated  with  the  total  strain  of  the  cyclic  stress  strain 

curve,  less  its  elastic  strain 

Ac i  ■*  Cl,,  “  Cj  —  2cr >/ £"j  (20) 


Estimating  the  propagation  strain  is  not  so  neat,  as  the  amount  of  propagation  depends  on 
the  value  of  GRF,  the  variable  being  solved  for.  Considering  the  term  0.3  m  e0  GRF,  it  is  not 
clear  what  to  use  for  c0  since,  in  effect,  e2  was  used  for  it  in  the  G2  term  of  GRF,  while  e,  was 
used  in  the  G 1  term.  Corresponding  values  of  this  differ  at  non-zero  stress  ratio,  as  detailed 
later.  To  obtain  a  closed-form  solution,  only  one  can  be  used.  In  this  case,  the  gradient  strain 
is  set  equal  to  the  maximum  strain,  corresponding  to  its  value  in  Gt.  This  is  the  more  accurate 
value  in  regions  where  the  term  dominates  the  G 2  term,  as  is  the  typical  of  higher  strength 
materials  of  high  environmental  sensitivity.  For  softer,  less  sensitive  materials,  the  propagation 


paths  are  dominated  by  the  G2  term  itself,  rather  than  G(M).  Here  this  inaccuracy  in  G(M)  is 
unimportant.  With  this  approximation,  then,  Eq.  (18)  becomes 


(21) 


(22) 
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or  for  the  parametric  M-family 


G(M) 


G:  +  (2*1  t,  G,  /0.15  m 


1  -  [2,w]  2  G] 

Similarly  to  Eq.  (17),  the  crack  growth  rate  for  this  family  is  simply 


da 

dN 


=  1.2  m  rT  G(M ) 


2.5  Corresponding  K-strain  Excursion 


(23) 


(24) 


Effects  of  load  ratio  R  on  the  values  of  G |  and  C2  are:  considered  next.  If  the  load  ratio  is 
zero  and  crack  tip  strain  is  proportional  to  the  stress  intensity  factor,  then  total  elastic  plus  plas¬ 
tic  strains  of  measurement  in  the  two  cycles  can  be  compared  directly.  However,  if  R  is  greater 
than  zero,  then  one  associates  monotonic  properties,  hence  G\,  with  a  strain  larger  than  that  for 
cyclic  properties,  hence  G2.  If  the  strain  excursion  is  taken  as  a  base,  as  is  AA\  then  the  basis 
of  association  is 

At  -  «2  -  (1  -  R)  t,  (25) 

This  means  that  to  calculate  a  GRF  for  a  given  value  of  strain  excursion  Ac,  one  uses  cyclic 

properties  at  strain  t2  =  At  to  compute  G 2,  while  using  monotonic  properties  at  strain  e,  = 

At/(1  —  R)  to  compute  the  corresponding  value  of  G\. 

What  level  of  stress  intensity  factor  excursion  AA"  is  required  to  effect  this  strain  At  at 
dT ?  The  model  allows  for  a  slight  degree  of  ineffectiveness,  or  slack,  in  coupling  the  overall 
K-excursion  to  the  strain  excursion  at  the  crack  tip.  This  slack  is  analogous  to  the  crack  closure 
concept  of  Elber  [7],  The  amount  of  this  slack  is  taken  as  one-half  the  monotonic  tensile  yield 
point  strain,  from  an  argument  viewing  the  compressive  stress-strain  quadrant  as  “off  limits"  for 
the  cyclic  strain  excursion  for  positive  stress  ratios  (4).  This  "closure  strain".  tCL  -  TYS/2E,  is 
taken  as  a  fixed  augmentation  of  the  monotonic  form  of  At  in  Eq.  (25).  In  the  cyclic  form,  it 
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is  added  only  to  the  extent  that  the  minimum  strain  R  Ae/(l-R)  is  less  than  the  closure  strain, 
or  the  difference  greater  than  zero.  Normally  this  term  vanishes  for  R  values  greater  than 


about  0.3.  With  this,  the  AX-proportional  strain  is  defined 

A  =  At  +  eCL  (26) 

“  *2  +  —  R  <2/(1  —  /»I>o  (27) 

*“  (1  —  R)  «i  +  tCL  (28) 

The  cyclic  excursion  AA  may  now  be  defined  after  Eq.  (6)  for  the  point  r  =  dT  =  2  rT 

AA  -  v/4^77  E  Ae*.  (29) 

An  illustration  of  foregoing  considerations  is  attempted  in  Fig.  1. 


2.6  Modification  of  Earlier  TLIM  Algorithm 

This  model  differs  from  that  of  Ref.  (4]  and  Ref.  [9]  in  one  major  respect,  as  well  as  in 
several  minor  ones.  The  major  change  is  that  the  gradient  strain  in  Gb  Eq.  (10)  and  in  G(M), 
Eq.  (21),  is  taken  as  a  maximum  strain  rather  than  a  strain  exco-ion.  When  the  strain  excur¬ 
sion  is  used  here,  the  level  of  G(A),  at  the  constand  N  fitting  a  stage  II  corrosion  fatigue 
growth,  rises  with  stress  ratio.  The  Bucci*  data  on  Ti-8-1-1  in  salt  water  appeared  to  exhibit 
this.  However,  when  this  is  applied  to  Vosikovsky  data  on  X70  steel  and  to  Ritchie  data  on 
300M,  the  R-effect  appears  erroneous.  It  is  now  believed  that  the  increasing  growth  rates  of 
the  Ti-8-1-1  data  should  be  associated  with  a  loading  in  which  A-™,,  exceeds  which  is 
Bucci’s  explanation  of  this  effect.  With  high  stress  ratio,  this  occurs  at  relatively  low  A K  levels, 
as  shown  later. 

The  numerous  minor  changes  are  to  be  noted.  In  Ref.  [4],  a  cyclic  stress  value  was  used 
in  calculating  as  well  G2.  The  gradient  strain  values  were  put  in  as  K -strains  e*  and  A«^, 


•References  to  authors'  works  noted  in  the  section  are  cited  later. 
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making  both  ordinate  and  abscissa  values  functions  of  the  closure  strain.  Finally,  the  constant 
in  G2  was  taken  as  unity  rather  than  VJ/2.  These  latter  two  changes  have  little  effect  except 
that  using  the  K-strains  seems  to  provide  slightly  more  stable  G(M)  mapping.  For  this  reason 
they  are  used  in  the  GRF  maps  to  be  shown  in  this  report.  However,  results  and  conclusions 
are  scarcely  affected  by  this  option.  For  future  work,  we  intend  to  use  the  simpler  algorithm  as 
put  forth  in  this  section,  and  so  recommend. 

One  further  variation  in  the  plots  to  be  shown  is  in  the  calculation  of  the  closure  strain. 
On  300M  steels  with  little  tempering,  the  normal  0.2%  offset  yield  point  is  markedly  higher 
than  the  proportional  limit.  To  reflect  this,  yet  not  destroy  the  usual  effectiveness  of  the 
TYS/2E  value,  the  closure  strain  was  entered  as  a  running  variable  eCL  =  tr^e^A/JE  Again 
the  effect  is  a  minor  one,  but  does  improve  correlation  in  a  few  cases. 

3.  PREPARING  OF  GRF  PARAMETRIC  CURVES 

To  compare  the  above-developed  algorithm  with  plotted  fatigue  crack  growth  data,  it  is 
necessary  to  construct  parametric  curve  sets  of  G(N)  and  G(M)  vs  Ae*.  This  task  is  too 
laborious  for  manual  calculation  and  graphing,  yet  it  is  readily  accomplished  with  a  suitable 
computer  program,  large  digital  computer  and  plotter.  Such  a  program,  available  from  the 
Naval  Research  Laboratory,  has  been  written  in  BASIC  for  the  NRL  DEC  PDP-10  computer. 
The  overall  procedure  involves  six  steps,  summarized  in  sections  to  follow. 

3.1  Measuring  Stress/Strain/Time  Curves 

The  TLIM  requires  a  mechanical  stress-strain  test  of  the  material  of  interest.  Specimens 
are  preferably  excised  from  used  fatigue  test  specimens,  with  longitudinal  axis  normal  the  frac¬ 
ture  plane;  ours  are  4.32  mm  (0.170  in.)  diameter,  with  about  9  mm  (0.35  in.)  length  of  uni¬ 
form  cylindrical  section  between  fillets.  End  buttons  1/2  -  20  NF  are  threaded  in  to  bottom 
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firmly  in  an  alignment  subpress,  then  subjected  to  tensile  and  compressive  deformation  in  the 
lower  head  space  of  10  KIP  (45  KN)  screw-driven  Instron  machine.  It  is  biased  with  heavy  coil 
springs  to  remove  backlash  in  load  reversal,  in  the  manner  of  Coffin  and  Tavernelli  [8], 

A  four-finger,  resistance-strain-gage-instrumented  clip  gage,  Fig.  2,  was  used  to  measure 
longitudinal  strain,  providing  a  stress-strain  curve.  It  is  designed  to  fit  outside  the  bi-lobed 
dimetral  strain  gage  used  in  previous  studies.  A  cross-plot  from  the  two  gages,  using  the 
second  pen  of  a  two  pen  X,  Yt,Y2  recorder,  provides  instant  measure  of  the  Poisson  ratio  (not 
shown  in  this  paper).  Fig.  3  shows  a  typical  stress-strain  curve  set. 

The  Instron  was  set  up  to  give  a  constant-load  strain  rate  of  about  10-3  sec.-1.  Data  for 
the  stress  relaxation  exponent  m  is  obtained  on  the  flat  region  of  a  cyclic  curve  by  arresting  the 
head  and  timing,  with  the  highest  available  amplification,  the  decay  in  load.  The  slope  of  a  log 
plot  of  load  vs  time  after  arrest  is  used  as  the  measure  of  m.  Usually  two  or  three  replicate 
tests  are  performed  on  a  given  material,  and  on  each  of  these  several  stress  relaxation  runs. 

3.2  Reading  and  Correcting  Flow  Data 

Values  of  stress  and  tangent  modulus  are  scaled  directly  on  the  stress-strain  curves  at 
selected  strain  stations.  The  curve-set  to  be  read  is  selected  from  those  available  after  visual 
comparisons  show  it  to  be  typical  and  the  record  of  good  quality.  Denoting  tensile  values  with 
subscript  T,  true  values,  used  in  the  equations,  with  none 


<T  t  ( 1  4  €  p) 

(30) 

0  j  4  O'  j  ( 1  4  €  y)  ( 1  -4“  C  7-)  ^ 

(31) 

In  (1  +  €  -p) 

(32) 

For  the  cyclic  curve,  stress  is  measured  from  the  compressive  toe,  converted  to  true  stress, 
then  divided  by  two  for  the  applicable  positive  portion  after  centroid  zeroing. 
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Typifying  records  used  in  (he  model,  this  shows  (he  full  monotonic  and  cyclic 
stress-strain  curves  for  2219-T851  Aluminum. 
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It  happens  that  the  strain-limited  growth-rate-factor,  G(M)  of  Eq.  (23),  is  extremely  sen¬ 
sitive  to  residual  elastic  strains;  they  should  be  completely  removed  by  use  of  Eq.  (20).  To 
assure  this,  the  measured  strain-station  value  is  replaced  by  one  calculated  from  measured 
stress  <r2  and  cyclic  Young’s  modulus  E2,  and  input  the  computer  program  as  true  values.  First 
cycle  properties  are  adequately  "trued"  in  the  program,  hence  may  be  input  as  measured 
engineering  properties. 


Monotonic  properties  of  structural  steel  usually  exhibit  an  upper/lower  yield  point  effect. 
The  ostensibly  zero  tangent  modulus  disrupts  the  model  by  producing  an  infinite  value  of  G,. 
On  the  other  hand,  good  correlations  are  obtained  by  removing  the  yield  point  effect.  Cer¬ 
tainly,  conditions  of  crack  tip  material  upon  cyclic  straining  do  not  allow  time  to  restore  the 
delayed  yield  effect,  as  by  strain  aging  processes.  The  measured  yield  point  effect  is  removed 
by  fitting  a  power  hardening  equation  to  the  region  beyond  that  of  the  lower  yield  point,  then 
extrapolating  with  it  back  under  the  measured  lower  yield  plateau.  The  set  of  equations  used 
follows: 


<t\  =  (t0  («,/«„)" 


0,  = 


1 


-1 


(33) 

(34) 


E 

<i  =  *p  +  vjE]  (35) 

The  strain  hardening  exponent  n  is  determined  as  the  slope  of  a  log  plot  of  (measured)  o,  vs  e, 


beyond  the  lower  yield  range,  and  E\  is  the  measured  monotonic  elastic  modulus. 


Although  zero  values  of  measured  0!  are  intractable,  the  measured  stress  values  produce 
no  infinity.  The  advantage  of  using  them  is  that  the  proportional  limit,  which  along  with  tcl 
sets  the  crack  growth  threshold  level,  can  be  determined  from  curve  measurement.  For  this 
paper,  values  of  monotonic  cr,  are  used  while  measured  values  of  9t  are  used  until  those  of  the 
power  law  estimate  exceed  the  measured  values.  This  occurs  at  a  strain  beyond  the  propor¬ 
tional  limit,  so  that  the  threshold  strain  can  be  identified  with  the  proportional  limit. 
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The  monotonic  and  cyclic  flow  data  is  listed  in  separate  data  file  "programs"  in  computer 
storage.  The  file  for  a  given  material  is  withdrawn  and  inserted  in  the  calculational  program 
using  the  BASIC-language  WEAVE  instruction  capability  in  preparing  the  program  to  RUN. 

3.3  Pairing  Monotonic  with  Cyclic  Properties  vs  Strain  Excursion 

Properties  from  the  two  curves  must  be  obtained  at  differing  strain  levels  «|  and  e2  for 
given  strain  excursion  Ac,  Eq.  (25).  The  stations  of  measurement  will  rarely  correspond.  To 
provide  correspondence,  a  log-uniform  set  of  strain  stations  is  established.  In  the  plots  to  be 
shown,  a  density  of  20  stations  per  decade  was  employed  to  approximate  the  density  of  meas¬ 
urement  points.  Greater  density,  requiring  more  computing  and  plotting  time,  is  of  little 
advantage.  Using  this  list,  we  interrogate  lists  of  cyclic  properties  vs  the  strain  e2  to  obtain 
these  properties  as  function  of  the  strain  excursion,  and  at  monotonic  properties  vs  the  strain  of 
€,(1-/?)  to  obtain  these  properties  as  function  of  the  maximum  strain.  Lists  of  values  are 
held  in  memory. 

3.4  Pairing  Monotonic  and  Cyclic  Properties  vs  A  K-strain 

The  fourth  step  augments  the  strain  excursion  values  with  the  closure  strain.  But  since 
the  augmentation  differs  for  each  set,  if  R  ^  0,  the  new  Ae*  strain  stations  are  no  longer  in 
registry.  Registry  is  restored  by  a  second  round  of  interpolations,  to  the  same  log-uniform 
strain-station  set,  now  stations  in  Ac*. 

It  should  be  noted  that  with  the  algorithm  described  here,  which  uses  a  fixed  value  of  the 
closure  strain  and  =  G(«|)  and  G2  **  G(t2),  it  is  possible  to  eliminate  the  second  round  of 
interpolations.  To  do  this  it  is  only  necessary  to  convert  measured  strain  values  «i  and  «2  to 
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At*  using  Eq.  (25-28).  However,  in  earlier  algorithms,  as  indeed  in  the  slight  variation  on  that 
described  above  which  was  used  to  generate  the  GRF  maps  of  this  report,  two  rounds  of  inter¬ 
polations  are  required. 

3.5  Plotting  Growth  Rate  Factors  Maps 

The  data  of  step  3.4,  or  of  step  3.3  if  only  one  round  of  interpolations  is  used,  is  then 
substituted  in  Eq.  (16)  for  G(N)  and  (23)  for  G(M).  On  the  G(N)  plot,  the  range  -21  ^  N 
<  -6  and  on  the  G(M)  plot,  -10  <  M  <  +6  is  found  to  suffice.  The  resulting  set  of  values 
of  G(N)  vs  A«*  and  G(M)  vs  Ae*  are  stored  in  "PLOT"  files  in  the  program,  which  are  called 
up  separately  for  execution  in  an  HP  7200A  Graphic  Plotter.  The  program  has  provisions  for 
selecting  the  scale  and  range  of  the  plot  to  match  the  data  plot,  as  well  as,  optionally,  drawing 
in  the  coordinate  scales,  as  shown  in  at  least  one  of  each  data  match  illustration  to  follow.  Here 
both  G(N )  and  G(M )  maps  have  been  superimposed  to  save  space.  The  computer  plots  of 
these  are  more  readily  distinguished  by  using  contrasting  pen  colors  for  each  one. 

3.6  Installing  Matching  Guidelines 

The  matching  process  involves  finding  what  values  of  the  fitting  parameters  A rr,  A rT/rT 
or  N,  Arr/rr/A«  or  M  provide  the  best  fit  of  the  data.  This  is  done  by  graphical  superposition 
of  GRF  vs  A e*  maps  upon  the  da/dN  vs  A K  data,  both  of  identical  (logarithmic)  scaling.  In 
doing  this,  it  is  convenient  to  install  match  guidelines  on  both  the  plots.  First  select  an  arbi¬ 
trary  point  on  the  crack  growth  data  plot:  da/dN  =  10  pm  cycle;  A K  =  10  MPa  Vm.  Solving 
Eqs.  (14)  and  (29)  using  these  coordinates,  plus  an  arbitrary  value  of  dT  —  100  ^m,  gives 

GRFref  -  1/6  m  (36) 

A «KR£F-  1000  (MPa)/ y/Trr  £, 

Other  values  of  dT  can  be  substituted  to  establish  a  line;  alternatively  one  of  slope  +  2  can  be 
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drawn  from  above-designated  reference  points,  by  regarding  the  2:1  dr-sensitivity  of  GRF  vs 
A«*.  The  matching  lines  shown  on  the  various  data-niatch  figures  follow  this  convention. 
Comparison  of  (7-maps  with  da/dN  data  is  by  superposition  of  the  respective  plots  with  match¬ 
ing  tines  held  coincident.  The  position  of  A t/cnet  on  the  AA-  scale  at  best  match  is,  con¬ 
veniently,  a  value  whose  square  is  the  size  of  the  dT  process  zone  in  microns. 

3.7  Number  of  Fitting  Parameters 

Regarding  the  number  of  fitting  parameters,  this  model  requires  from  one  to  five,  these  of 
three  different  kinds.  The  size  parameter  dT  is  always  required  If  there  is  little  or  no  environ¬ 
mental  effect,  the  data  should  correspond  to  G(N)  «=  (7ronly,  for  a  one-parameter  (dr-only) 
fit.  If  sensitivity  to  the  air  environment  exists,  data  will  generally  correspond  to  a  curve  of  con¬ 
stant  M  =  M,n  for  a  two-parameter  fit.  In  aggressive  liquid  environments,  two  additional 
stages  must  be  characterized,  stage  1  growth  by  another  higher  value  of  M  «  A/;;  and  stage  II 
bridging  I  to  III  with  a  path  of  constant  N  =  N „  for  a  four  parameter  fit.  Finally,  at  high  stress 
ratio,  one  can  detect  a  "knee"  protruding  down  from  the  constant  —  M/  stage  I  trace  which 
appears  to  be  matched  by  a  curve  of  very  low  N  «=  NTH  value.  The  corresponding  depth  of 
environmental  attack  is  miniscule,  hence  not  unreasonable.  But  with  it,  a  maximum  total  of 
five  fitting  parameters  are  needed.  However,  the  same  parameters  apply  to  all  positive  stress 
ratios  without  change. 

The  two  curve-families  differ  in  frequency  sensitivity.  Curves  of  constant  M,  like  the 
stages  I  and  II  of  CFCG  which  they  track,  tend  to  be  frequency-independent,  presumably 
because  here  the  amount  of  surface  attack  is  limited  by  the  degree  of  surface  distortion,  vis  a 
vis  plastic  deformation  of  the  ligament.  On  the  other  hand,  the  curves  of  constant  i V„,  as  the 
stage  II  they  track,  tend  to  be  frequency/cycle-duration  dependent,  presumably  related  to  the 
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time  dependency  of  surface  attack  or  corrosion.  One  may  thus  regard  growth  in  stages  1  and  III 
as  of  strain-limited  surface  attack,  rather  than  time-limited,  and  stage  II  as  of  time-limited  sur¬ 
face  attack,  rather  than  strain-  limited. 

4.  CASE  STUDIES 


The  substantial  number  of  case  studies  now  to  be  displayed  is  possible  due  the  generosity 


of  authors  who  donated  specimen  material  and  permission  to  attempt  this  modeling.  Table  I 
provides  a  listing  of  alloy  compositions.  Table  II  provides  a  list  of  heat  treatments  and  resulting 
grain  size  and  mechanical  properties,  and  reference  to  the  FCGR  data  analysis  figures  of  this 


text.  Each  of  these  figures  contains  the  two  families  of  growth  rate  factor  maps,  superimposed, 


from  which  the  fitted  curves  were  selected,  and  then  the  FCGR  data  with  selected  curves  over¬ 
laid.  The  characterizing  values  of  dT,  NTH,  Mh  N„  and  Mllh  as  appropriate,  are  shown  on  each 
figure,  but  not  elsewhere  tabulated  (except  for  dT).  However  some  values  of  surface  attack  per 


I.  Table  of  Alloy  Composition 


Alloy 

Base 
+  wt% 

c 

Si 

— 

Mn 

— 

Ni 

Cr 

Mo 

V 

Co 

Cu 

Al 

Nb 

S 

P 

A  36 

Fe 

0.17 

0.29 

0.86 

0.03 

0.021 

0  015 

BS  4360- SOD 

Fe 

0.17 

0.358 

1.35 

0.07 

0.09 

0.01 

0.17 

0.05 

0.03 

0023 

0.037 

X  65 

Fe 

0.16 

0.33 

1.34 

0.046 

0.031 

0  009 

0.006 

X  70 

Fe 

0.06 

0.28 

1.90 

0.25 

0.39 

0.073 

0.058 

0.009 

0.002 

SM  580 

Fe 

0.14 

0  32 

1.30 

0.024 

0.020 

0023 

0.038 

0.006 

0,022 

HT  80 

Fe 

0  13 

0  19 

0.94 

0.025 

0.46 

0.14 

0.027 

0.051 

0.007 

0.014 

HY  130 

Fe 

009 

0.26 

069 

5.0 

0.48 

0.47 

0.08 

0.05 

0.05 

0.007 

0009 

10  Ni 

Fe 

0  12 

0.07 

0.28 

10.29 

7.03 

1.03 

8.07 

0.006 

0.008 

4340 

Fe 

0.40 

0.31 

0.72 

1  85 

0.82 

0.22 

0.020 

0.010 

300M 

L  Fe 

042 

1.59 

0.76 

1.76 

0.76 

0.41 

0.10 

0.002 

0.007 

Base 

r- 

+  wt% 

Al 

V 

Me 

Fe 

C 

O 

N 

H 

R  14A 

Ti 

6.0 

4.1 

0.05 

0.023 

0.06 

0.008 

00050 

R  55A 

Ti 

6.1 

4.0 

0.18 

0.02 

0.11 

0.009 

0.0069 

R  54A 

Ti 

6.6 

4.4 

0.20 

0.02 

0.18 

0.014 

0.0058 

R  23C 

Ti 

6.7 

4.3 

0.10 

0.03 

0.20 

0.011 

0.0060 

W&D  6-4 

Ti 

6.3 

3.9 

0.07 

0.01 

0.14 

0.01 

0.005 

R7 

Ti 

7.4 

0.86 

1.16 

0.14 

0.08 

0.09 

0.004 

R  64  A 

Ti 

7.8 

II 

10 

0.07 

0.03 

0.11 

0.015 

0.0046 

WAD  8-1-1 

Ti 

7.7 

1.0 

1.0 

0.08 

0.022 

0.09 

0.009 

0.005 

Base 

+  wl% 

Cu 

Zn 

Mg 

Si 

Mn 

Ti 

Cr 

Fe 

2124 

Al 

3.9 

1.3 

0.08 

0.54 

<0.03 

0.14 

2219 

Al 

6.28 

0.025 

0003 

0.088 

0.25 

0.051 

<0.0001 

0.25 
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il.  Table  of  Alloy  Heat  Treatment,  Grain  Size,  Mechanical  Properties 


Fig- 

Alloy 

l  bKA 

Dir 

<GPa> 

<GP*J 

TYS°* 

KJPat 

UTSU’ 

•» 

m 

Gram 

Steel 

t 

t~  :  ,-T-  :=z=,,J 

•  ' 

Si/e 

...  : 

4 

A?6 

4R 

TL 

207 

0  26/ 

0  48/ 

18 

0  0072 

77 

5 

HS4160  5l)l> 

18 

AK 

:  207 

0  37/0  37 

0  53/0  54 

0  26 

0  0078 

87  f 

6 

X65 

t) 

API  Spec  51. \ 

TL 

207 

0  45/0  46 

0  59/0  57 

=0  2 

0  0064 

18 

25 

7 

X?0 

]9 

C  rinif.l-Ri.lled.  API  Spec  5LX 

LT 

207 

0  46/0  53 

0  63/0  67 

0  16 

0  0065 

4  5 

20 

8 

sm  *8  y 

12 

V3U\  wq.  6M»  I  AC  1 

LT 

i  210 

0  59/0  59 

0  68/0  67 

0  09 

0  0059 

18 

9 

HT  80 

12 
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cycle,  as  calculated  from  the  figures,  are  plotted  in  summary  figures.  The  figure  captions  con¬ 
tain  reference  to  the  author  of  the  FCGR  data,  and/or  to  the  donor  of  specimen  material.  The 
figures  are  ordered  in  the  grouping  and  sequence  of  Table  II.  Some  notes  on  particular  cases  of 
the  data  collection  follow. 

4.1.  Iron  Based  Alloys 

Figures  4  to  16  display  TLIM  analysis  of  the  sixteen  steels  of  differing  composition  and/or 
heat  treatment. 

Linear  elastic  fracture  mechanics  parameters  can  be  applied  to  mild  steels  because  fatigue 
crack  growth  occurs  at  A  A-  levels  well  below  Kk.  Actually  the  fatigue  instability  limit,  the 
AK  for  indefinitely  rapid  FCG  rate,  is  generally  much  less  than  K,c.  In  the  TLIM,  this  is  asso¬ 
ciated  with  cyclic  hardening,  which  reduces,  by  as  much  as  a  factor  of  2  to  3,  the  cyclic  strain 
for  maximum  (tensile)  load.  Thus  even  in  a  steel  as  soft  as  A-36,  Fig.  4,  a  substantial  part  of 
the  growth  curve  can  be  measured,  and  fitted  by  TLIM  procedures.  However  at  a  stress  ratio  as 
high  as  0.8,  only  the  threshold  region  can  be  interpreted  with  K-parameters.  Here  the  general 
shape  of  the  threshold  plateau  does  appear  to  be  roughly  followed  by  the  NTH  =  — 18  curve.  A 
more  refined  closure  strain  algorithm  is  needed  for  more  accurate  prediction  in  this  region. 
The  manner  of  removing  upper  yield  point  effects  is  a  factor  here. 

The  TLIM  match  of  BS4360  Grade  50D  in  both  room  air  and  5-10  degrees  Centrigrade 
sea  water,  Fig.  5,  seems  rather  satisfactory,  as  also  that  for  line  pipe  steels  X65  and  X70  in  air 
and  room  temperature  salt  water.  Figs  6  and  7.  In  X65,  Fig.  6,  the  shift  to  the  right  with  fre¬ 
quency  of  the  low  H2S  crude  data  could  be  an  effect  of  strain  rate  on  the  size  and  shape  of  the 
stress-strain  curve,  to  be  discussed  later. 
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The  extensive  R-effect  data  in  air  on  SM58Q  and  HT80  steels,  Figs.  8  and  9  is  reasonably 
well  matched.  The  TLIM  curves  for  R  —  0  does  separate  from  the  high-R  curves,  in  mid¬ 
range  A K,  more  than  the  data,  particularly  in  the  SMS8Q.  Here  the  modification  of  the  TLIM 
equations  to  include  the  closure  strain  in  the  gradient  strain  term  of  G,  and  G2  tends  to  exag¬ 
gerate  this  discrepancy. 

Regarding  HY130  and  the  10  Ni  steel,  Figs.  10  and  11,  there  is  some  difficulty  in  matching 
both  ends  of  the  crack  growth  range  with  the  growth  rate  factor  maps.  If  one  fits  to  favor  the 
threshold  region,  as  has  been  done  here,  the  run-out  in  the  fatigue  instability  region  falls  short; 
and  vice  versa.  The  behavior  is  much  as  though  the  closure  strain  term  is  too  large.  It  is  possi¬ 
ble  that  residual  stress  from  the  quenching  operation  in  heat  treatment  could  have  an  effect. 
Significant  influences  of  this  kind  have  been  reported  by  Sullivan  and  Crooker  [29]  and  by 
Ritchie,  Smith  and  Knott  [30].  Another  problem  is  that  the  monotonic  tensile  stress-strain 
curves  of  these  materials  is  so  extraordinarily  flat  that  it  is  difficult  to  determine  the  point  of 
maximum  load.  Large  effects  of  strain  rate,  via  cyclic  frequencies,  are  likely  here.  In  view  of 
these  uncertainties,  present  matching  parameter  values  on  these  two  steels  should  be  regarded 
as  provisional. 

If  residual  stress  is  a  factor  in  the  high  strength  quenched  and  tempered  steels,  what  about 
4340  and  300M?  In  the  4340,  Fig.  12,  the  data  is  not  critical  as  the  growth  threshold  region  has 
not  been  characterized,  hence  the  model  is  easily  fitted  to  mid  and  upper  ranges.  In  the  300  M 
a  full  range  of  FCG  data  is  available.  Indeed  the  higher  tempering  temperatures,  hence  lower 
residual  stress,  give  the  better  fits,  particularly  in  the  threshold  region.  A  point  of  caution  here 
is  that  the  tensile  specimens  of  300M  were  heat  treated  separately  on  plate  material  provided  by 
Ritchie;  they  were  not  cut  out  of  heat  treated  fracture  specimens,  as  is  our  preferred  practice. 
It  appears  that  this  has  failed  to  develop  the  strength  properties  of  the  test  specimens,  as  seen 
in  Table  II,  where  the  two  values  are  compared. 
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4.2.  Titanium  Based  Alloys 

Figures  17  to  27  display  TLIM  analysis  of  the  seven  Ti-6A1-4V  alloys  and  four  Ti-8A1- 
lMo-lV  alloys  of  differing  composition  and/or  heat  treatment. 

The  Yoder  data  on  Ti-6-4,  as  on  Ti-8-1-1,  was  meant  to  characterize  changes  in  FCGR  as 
a  function  of  microstructural  and  interstitial  concentration  changes.  Hence  the  amount  of  data 
collected,  though  carefully  measured,  was  only  as  extensive  as  required  to  characterize  the 
mid-range  growth  rates.  This  includes  the  mid-range  slope  inflection,  which  Yoder  associates 
with  correspondence  of  grain  size  and  the  reverse  plastic  zone  size.  With  this,  the  TLIM  curve 
fit  is  not  as  rigorously  bounded  as  with  more  extensive  data  sets.  However  the  curve  matching 
seems  fairly  tight,  and  some  sort  of  inflection  point  is  suggested  directly  from  the  shape  of  the 
constant- Mw  curves. 

The  two  data  collections  from  Wanhill  and  Doker,  Figs.  23  and  27,  touch  on  the  question 
of  two  stage  hardening  effects.  With  continued  cycling,  the  strain  hardening  rate  of  many 
titanium  alloys  tends  to  decrease,  then  at  about  half  of  the  strain  excursion,  increase  to  recover 
the  earlier  strength  loss.  It  was  thought  that  the  abrupt  jump  in  growth  rate  shown  by  these 
authors  might  be  associated  with  this  effect.  To  this  end,  the  second  full  cyclic  excursion, 
showing  a  moderate  but  not  fully  developed  two  stage  hardening,  was  measured.  The  second 
cycle  map  does  appear  to  provide  a  closer  fit  of  the  FCGR  data.  However,  first  cyclic  curves 
also  fit  reasonably  well  and  since  these  had  been  used  in  all  other  cases,  the  dT  value  from  its 
match  is  used  in  comparisons  of  fitting  parameters  to  follow.  This  is  only  a  token  "glimpse"  of 
the  subject,  much  more  should  be  done  on  it. 
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matched  to  air  FCGR  data  of  Yoder,  Cooley,  and  Crooker.  Ref.  18  and  19 
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Finally  in  the  titanium  alloys,  results  on  the  R-7  alloy  used  by  Bucci  is  shown  in  Fig.  24. 
One  need  no  longer  be  concerned  that  the  growth  rate  in  dry  argon  requires  an  M  value  com¬ 
parable  to  those  for  room  air  environment.  Later  work  of  Wei  and  Ritter  (31)  show  that 
minute  traces  of  water  could  have  caused  this.  The  salt  water  data  is  nicely  matched  by  /V/,  = 
-7.5  up  to  the  region  in  which  Km.i%  exceeds  whereupon  a  marked  acceleration  is 

observed. 

4.3.  Aluminum  Base  Alloys 

Figures  28  and  29  display  TL1M  analysis  of  the  two  aluminum  alloys. 

Regarding  the  2I24-T851  data.  Fig.  28,  the  mid-  to  upper-range  fit  is  fairly  satisfactory. 
However  the  TLIM  estimates  fall  well  below  measured  rates  in  the  threshold  region.  It  is  here, 
much  as  with  the  high  strength  steels,  as  though  the  closure  strain  is  too  large.  In  the  2219- 
T851,  Fig.  29,  the  fit  is  better.  In  both  alloys  thin  specimens  were  cut  out  of  thick  plates  so  it 
is  unlikely  that  surface  stresses  due  to  quenching  from  solution  anneal  played  a  role  in  either 
case. 

In  the  2219  data  plot.  Fig.  29,  an  attempt  has  been  made  to  model  the  data  at  stress  ratio  -I. 
It  is  argued  that  one  effect  of  compressing  a  crack  might  be  to  overcome  the  resistance  of  crack 
closure.  To  simulate  this  in  the  TLIM,  one  can  remove  the  closure  strain  term  and  run  the  cal- 
culational  program  with  R  set  to  zero.  The  resulting  curve  seems  to  fit  the  data  rather  nicely. 
In  fact  in  the  high  A K  region,  the  fit  is  better  than  for  the  R  =  0.1  data,  where  the  runout 
rates  exceed  the  TLIM  curve,  and  oddly,  the  R  =  —  1  data  trend. 

In  this  data,  only  the  positive  part  of  the  stress  excursion  is  used  to  calculate  AA\  even 
though  the  minimum  stress  is  negative.  When  the  Sasaki  et  al  data  (13)  is  treated  in  this  way 
(not  shown  in  Figs.  7  and  8),  the  TLIM  estimate  is  too  far  to  the  left,  as  though  some  of  the 
closure  strain  term  should  have  been  retained.  More  work  is  needed  on  this  subject. 
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Fig.  28  —  GRF  curves  for  2124-T851  aluminum  at  iwo 
stress  ratios  is  matched,  favoring  the  higher  levels  of 
A K,  to  air  FCGR  data  of  Ruschau,  Ref  25. 
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5.  COMPARISON  OF  FITTING  PARAMETERS 

It  is  suggested  that  the  degree  of  correspondence  between  TLIM  contours  of  fixed  param¬ 
eter  designations  and  the  fatigue  crack  growth  rate  data.  Figs.  4-29,  is  encouraging.  Some  fits 
are  extremely  close;  a  few  leave  something  to  be  desired,  though  not  without  extenuating  cir¬ 
cumstances.  Since  most  of  the  correlations  are  close,  it  could  be  informative  to  consider  impli¬ 
cations  of  trends  in  fitting  parameter  values. 

5.1,  Trends  in  Process  Zone  Size 

The  process  zone  size  is  basic  to  the  TLIM.  It  is  required  as  a  fitting  parameter  for  every 
sort  of  environmental  condition,  including  the  dry  inert  gas  or  high  vacuum.  Indeed  some  sort 
of  size  parameter  is  generally  necessary  in  any  model  relating  flow  properties  to  fracture 
behaviors.  Recall  here,  that  the  process  zone  size  is  the  distance  ahead  of  the  crack  tip  in  the 
inverse  square  root  strain  singularity  that  provides  the  best-fit  proportionality  factor  between 
measured  strain  in  the  tensile  test,  plus  the  closure  strain,  and  AA  or  A/c  in  the  fracturing  test. 
The  value,  shown  on  each  data  match  figure,  is  obtained  from  the  intercept  of  the  reference 
point  XeKRFF  on  the  AA  scale.  The  square  of  this  value  of  AA  in  MPa  is  dT  in  juin.,  as  noted 
earlier. 

Yoder  has  investigated  effects  of  microstructural  modification  on  FCGR  in  several 
Titanium  alloys.  He  finds  a  substantial  lowering  of  growth  rate  in  mid-range  with  increase  in 
effective  grain  size.  In  beta-annealed  microstructures,  grain  size  is  identified  with  the  Widman¬ 
statten  packet  size  whereas  in  lower  temperature  anneals  it  is  associated  with  the  alpha  phase 
grain  size.  Heat  treatments  and  average  grain  size  values  are  given  in  Table  II,  and  data 
matches  leading  to  dT  in  Figs.  17-23,  25  and  26.  Fig.  30  is  a  (log)  cross  plot  of  the  dT  process 
zone  size  against  mean  grain  size  of  the  titanium  alloys.  It  shows  a  reasonable  proportionality 
between  the  variables,  perhaps  a  better  correspondence  with  the  process  zone  radius  than  its 
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diameter.  This  correlation  suggests  that  one  might  associate  the  improvement  in  FCGR  resis¬ 
tance  due  to  a  coarse  microstructure  to  increased  process  zone  size.  Roughly  speaking,  if 
environmental  factors  are  similar,  increases  in  dT  transfer  G  values  of  a  fitting  TL1M  curve 
along  the  match  line  direction,  in  ratio  2:1;  e.g.,  a  20%  increase  in  d,  increases  growth  rate  by 
20%  percent  but  the  A  A  to  effect  it  is  increased  by  10%.  If  the  slope  of  the  log  daldN  vs  A  A 
curve  is  greater  than  2,  as  true  of  these  alloys,  the  net  effect  is  a  decrease  in  growth  rate  for 
given  AA.  Some  further  decrease  could  result  from  changes  in  strength  level  as  affecting  the 
form  of  the  stress  strain  curves.  As  noted  earlier,  the  value  of  dr  for  the  Wanhill  and  Doker 
Ti-6A1-4V  is  that  for  a  first  full  cyclic  stress-strain  curve,  as  is  the  case  for  all  others,  not  the 
second  cyclic  curve  which  is  matched  in  Fig.  23. 

The  correspondence  of  dT  with  grain  size  suggests  the  plotting  of  it  against  yield  strength, 
to  which  grain  size  is  often  related.  Cross  plots  vs  the  0.2%  offset  monotonic  yield  strength  are 
shown  in  Fig.  31  for  the  steels  and  in  Fig.  32  for  the  titanium  alloys,  in  the  uppermost  data 
band.  The  process  zone  size  in  steel  appears  in  steady  decline  with  increasing  yield  strength. 
The  slope  of  the  band  overlaid  is  as  though  a  Hall-Petch  relationship  governed;  i.e.,  </,,  via 
grain  size,  proportional  to  minus  two  power  of  yield  strength,  or  a  log  slope  of  -2. 

In  the  titanium  alloys,  the  overlaid  steel  trend  band  from  Fig.  31  seems  to  bound  the  data. 
However  here,  the  trends  do  not  suggest  a  -2  power  slope.  In  the  6-4  data,  increasing  oxygen 
content  increases  yield  strength  with  little  effect  on  the  d, -value,  as  one  might  expect  of  an 
interstitial  hardener.  But  the  increase  in  yield  strength  with  reduced  d,  (hence  grain  size'.’)  is 
far  smaller  than  seen  in  the  steels.  Conversely,  the  loss  of  strength  with  improvements  due 
increased  dr,  is  far  less  than  might  be  feared.  This  lower  sensitivity  is  also  seen  in  the  data  on 
the  8-1-1  alloy  on  Fig.  32. 
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Eig  3 2  —  The  si/e  ami  surface-attack -depth  values  derived  from  T1  IM  matching  ot  the 
various  titanium  base  alloys. 
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5.2.  Trends  in  Frequency-Independent  Parameters 

The  normal  frequency-independent  parameters  of  the  TLIM  are  two:  the  values  of  stage 
111  strain-limited  surface-attack  Mm  corresponding  to  air  environment  FCG,  and  characterizing 
stage-I  aqueous  environment  CFCG  A//.  In  the  sea-level  air  environment,  adequate  time  for 
saturation  of  the  surface  attack  process  is  thought,  after  results  of  Bradshaw  and  Wheeler  [26], 
to  be  far  shorter  than  any  fatigue  cycle.  In  the  aqueous  environments,  the  upper  limit  of  stage 
I  growth  would  normally  be  a  function  of  cyclic  frequency.  It  appears  also  that  the  threshold 
level  of  surface  attack,  NTH  or  (A rT)TH  also  may  be  frequency  independent,  unlike  stage  II 
attack  in  corrosion  fatigue,  which  is  frequency  dependent,  as  illustrated  later. 

Recall  that  values  of  these  fitting  parameters  are  obtained  from  the  data  match  plots, 
where  two  to  the  N  or  M  is  the  value  of  [A rT/rT]  or  [A/>/rr/A«]  respectively  for  the  fit.  Since 
the  value  of  rT  is  determined  by  fit  location,  the  absolute  value  of  A rT  or  Arr/Ae  is  simply  the 
product  or  2  U  times  rT.  Such  values  as  could  be  obtained  are  plotted  vs  yield  strength  for  the 
air  data  in  Figs.  31  and  32. 

Values  of  (Arr/A«)w  are  too  small  for  determination  in  the  softer  steels.  In  these,  the 
air  data  is  essentially  a  (?2-only  match;  negligible  air-environment  effect.  However,  in  the 
higher  strength  4340  type  steels.  Figs.  12-16,  measurable  values  are  discernible  —  of  the  order 
of  1/2  micron  per  unit  strain. 

The  salt  water  environment,  particularly  at  higher  levels  of  impressed  potential,  does  give 
measurable  values  of  stage  I  strain-limited  surface  attack  (Arr/A«),  in  the  mild  steels.  Values 
are  plotted  in  Fig.  31  for  the  extreme  cases:  free  corrosion  in  -0.65V  and  closely  coupled  zink 
potential,  —  -  1.1V.  These  lie  above  the  band  laid  over  the  air  values,  as  expected.  Unexpect¬ 
edly,  they  indicate  a  downward  trend  with  increasing  yield  strength,  almost  corresponding  to 
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that  for  dT.  It  is  almost  as  though  the  relative  surface  attack  \rT/rr/ A«  is  the  more  invariant 
parameter.  The  trend  in  (\rr/kt)  m  air  values  also  suggest  this  trend,  but  not  so  distinctly. 

The  residual  attack  (A rT)TH,  the  lowest  band  in  Fig.  31,  appears  only  the  order  of  a  lattice 
spacing  in  the  softer  steels,  rising  some  factor  of  ten  for  the  hardest  steels. 

Regarding  the  titanium  alloys.  Fig.  31,  stage  111  surface  attack  (\rr/be)„,  appears  of 
similar  severity  to  that  for  the  steels  but  with  the  Ti-6-4  alloy  definitely  less  susceptible  than  the 
Ti-8-1-1.  No  threshold  data  was  available  to  permit  estimate  of  (A  rT)TH  for  these  alloys. 

5.3.  Trends  in  Frequency-Dependent  Parameters 

The  frequency-dependent  parameter  of  the  TL1M  is  the  value  of  N,  or  corresponding  A rr 
per  cycle,  fitting  stage  II  corrosion  fatigue  growth  region.  Such  values  as  are  available  in  this 
collection  are  shown  in  Fig.  33.  Here,  on  a  logarithmic  plot,  the  value  of  A rr  is  plotted  against 
the  loading  time  of  the  cycle.  Using  loading  time  rather  than  the  half  or  full  period  follows  the 
finding  of  Barsom  [27]  and  of  Kawai  and  Koibucci  [28]  that  this  is  the  governing  time  in  corro¬ 
sion  fatigue  below  Kht.c.  Loading  time  is  taken  as  0.35  of  cycle  period  for  sinusoidal  wave  form 
and  0.5  for  triangular,  the  latter  form  due  to  Vosikovsky.  The  slightly-tempered  4340  steels,  M 
and  N  in  Fig.  33,  are  also  plotted  this  way  even  though  the  cycling  is  at  a  Amay  generally- 
exceeding  the  A 'ivc  threshold,  for  which  a  longer  period  fracf'  ,n  would  be  appropriate.  In  the 
case  of  the  R-7  titanium  of  Fig.  24,  the  full-cycle  period  would  seem  warranted  for  the  regions 
in  which  Kmm  ^  K,ur. 

In  evaluating  a  plot  as  Fig.  33,  one  should  place  less  confidence  in  the  low  values  of  A rt, 
as  the  small  differences  between  M,  and  Mm  paths  in  this  region  make  it  difficult  to  lit  Also 
with  small  values  of  the  environmental  term,  we  should  be  in  the  transition  region  between 
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cyclic  and  monotonic  curve  shape  control,  and  the  crudity  of  the  abrupt  transition  algorithm 
could  lead  to  inaccuracy.  A  slope  of  plus  1  on  this  plot  would  suggest  a  constant  value  of  sur¬ 
face  attack  rate,  called  in  earlier  papers.  There  appears  to  be  a  tendency  for  Vt  to  decrease 
with  cycle  period  at  longer  cycle  dura'  ons,  presumably  due  to  some  sort  of  corrosion  saturation 
effect. 

6.  DISCUSSION 

A  few  caveats  and  generalizations  are  put  forth  in  what  follows. 

6.1.  Sources  of  Correlation  Errors 

There  are,  of  course,  at  least  two  general  sources:  those  of  the  modeling  and  those  of  the 
crack  propagation  data.  The  model  is  admittedly  a  primitive  one.  Its  basic  hypothesis, 
constant-load-at-dr  maintenance,  can  only  be  viewed  as  a  possible  necessary,  but  surely 
insufficient,  condition.  Properties  of  uniformly  strained  uniaxial  specimens  are  applied  to  a 
region  of  high  strain  gradient  and  triaxial  stress  state.  The  characterizing  tensile  specimens  are 
enormous  compared  to  the  to  the  zone  of  interest,  which  is  in  fact  of  the  size  scale  of  some  of 
the  strength  affecting  microconstituents.  Indeed,  if  inclusions  are  stringerlike,  this  far  exceeds 
the  process  zone  size.  The  division  of  roles  of  monoionic  vs  cyclic  properties  is  abrupt;  the 
crack  closure  term  rudimentary. 

There  are,  additionally,  errors  of  recording  and  measuring  the  stress-strain  curves  which 
can  affect  accuracy.  The  threshold  determination  is  extremely  sensitive  to  curve-shape  around 
the  proportional  limit.  The  estimated  points  of  fast  fracture  instability,  or  fatigue  instability,  is 
sensitive  to  curve  shapes  around  maximum  load.  Aside  from  instrument  errors,  these  things 
are  affected  by  strain-rate  and  cyclic  hardening/ softening  effects.  The  "two  stage"  hardening 
tendency  in  cycled  titanium  alloys  is  of  concern. 
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Regarding  errors  in  growth  rate  data,  as  sophistication  in  FCG  measurements  grows,  more 
sources  of  data  variability  are  being  discovered.  False  thresholds  pervade  the  literature  — 
hopefully  none  in  the  data  used  here.  Work  by  Crooker  and  Sullivan  [29]  show  large  effects  of 
residual  stress  in  the  specimen.  Actually,  standard  ASTM  methods  for  FCG  measurement  are 
only  now  being  issued.  Thus  it  is  not  alone  the  model  imperfections  which  must  underlie  all 
discrepancies. 

6.2.  dT  vs  Microstructural  Sizes 

The  1962  work  of  Birkle,  Wei  and  Pellissier  [32]  associated  the  dr  process  zone  size  with  the 
average  spacing  of  largest  inclusions  of  manganese  sulfide  in  4345  steel.  Thus  it  is  somewhat 
unexpected  that  the  present  results  on  titanium  alloys,  and  by  strong  inference,  those  on  struc¬ 
tural  steels,  indicate  an  association  with  average  grain  size.  In  the  titanium  alloys,  inclusions 
could  hardly  be  a  factor,  as  few,  if  any,  are  present.  It  is  possible  that  in  very  fine  grained 
materials,  such  as  quenched  and  tempered  4345,  control  of  the  process  zone  size  shifts  from 
grain  size  to  a  coarser  microconstituent  such  as  the  inclusions.  Further  work  is  needed  here. 

6.3.  dT  from  Low  Temperature  Klr(k)  Matchings 

Another  conflict  with  earlier  concepts  is  in  the  magnitude  of  dT  values.  Typically  here  for 
structural  steels,  dr  values  from  10  to  100  microns  are  required  for  matching.  Yet  in  extensive 
low  temperature  matching  of  the  isothermal  strain  hardening  exponent,  or  strain  for  maximum 
load,  with  plane  strain  fracture  toughness,  values  on  order  of  magnitude  smaller  were  found 
[33],  It  is  clear  that  the  latter  values  are  too  small  to  explain  the  room  temperature  fracture 
toughness  of  steels.  One  possible  explanation  is  that  low  temperature  toughness  is  to  be  associ¬ 
ated  with  the  adiabatic  strain  hardening  rate,  much  lower  than  isothermal  at  the  reduced  specific 
heats  of  cyrogenic  temperatures.  Here  although  the  process  zone  would  be  isothermally 
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quenched  for  any  feasible  crack  loading  time,  the  existance  of  an  instability  path  to  trigger  a 
lower  energy  cleavage  propagation  leaves  a  choice.  Given  the  general  perversity  of  nature,  and 
the  unstable  nature  of  cra  'k  propagation,  the  option  for  the  easier  path  for  fast  fracture  initia¬ 
tion  is  within  the  realm  of  possibility. 

This  concept  poses  a  problem  at  temperatures/loading  rates  for  which  shear  rather  than 
cleavage  processes  characterize  the  separation.  For  example,  in  some  aluminum  alloys,  fracture 
toughness  actually  rises  at  cryogenic  temperatures  [34]  as  does  the  isothermal  strain  hardening 
exponent,  yet  the  decline  of  specific  heat  is  not  atypical  for  this  metal. 

An  alternate  to  the  constant  dT  concept  is,  after  Hahn  etal  [35],  is  a  process  zone  limited 
by  intercleavage  facet  ligamentary  material  to  a  size  much  smaller  than  that  of  the  whole  grain. 
Then  in  the  transition  to  shear-dominated  fracture,  some  sort  of  a  transition  in  microseparation 
mode  occurs,  shifting  the  process  zone  size  to  control  by  grain  size.  Spitzig  [36]  supports  an 
explanation  of  this  sort. 

6.4.  Cleavage  in  FCGR  Prediction 

One  should  be  cautious  in  attempting  to  use  the  TLIM  for  FCGR  predictions  when 
cleavage  is  included  in  the  fracture  process.  This  model  implicitly  assumes  a  ductile  tear  type 
instability  as  a  local  process.  The  additional  growth  due  to  local  bursts  of  cleavage  are  seen  by 
Stonesifer  [37]  and  Richards  and  Lindley  [38]  to  accelerate  the  average  crack  propagation  rate. 
The  TLIM,  as  presently  constituted,  has  no  way  of  accounting  for  this. 

Stonesifer  finds  the  fatigue  upper  limit  Kf(  to  drop  toward  the  threshold  lower  limit  as 
temperature  is  reduced.  He  associates  this  with  the  drop  in  K/c  with  decreased  temperature. 
Since  this  is  cleavage  fracture,  the  level  of  Klc  cannot  be  predicted  from  non-cleavage  values  of 
dT,  as  discussed  in  section  6.2,  except  possibly  as  an  adiabatic  initiation  trigger. 
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While  cleavage  fracture  is  not  handled  here,  the  model  does  seem  to  provide  useful  corre¬ 
lations  in  regimes  of  FCG  characterized  by  a  variety  of  other  non-cleavage  fractographic 
features.  For  example,  near  the  growth  threshold,  fracture  surfaces  can  be  characterized  by 
almost  planar  separation,  sometimes  modulated  by  waves  of  minute  striations.  Scant  evidence 
here  of  the  dT  ligament  such  as  the  ductile  tear  dimple. 

6.5.  Yield  Strength  Effect  on  KTH;  on  KIc 

The  model  algorithm  permits  a  direct  estimate  of  the  threshold  for  FCG.  In  the  non¬ 
vacuum  environment,  it  is  proportional  to  the  monotonic  yield  point  strain,  from  Eq.  (26)  and 
(29)  approximately 

A Kth  =  (1.5  -  R)  TYS  y/2^d~T  (37) 

This  is  a  slight  overestimate,  as  the  program  detects  the  proportional  limit,  rather  than  0.2% 

yield  stress.  Present  results  for  steels.  Fig.  3 1 ,  show  a  decline  in  dT  with  yield  strength,  approx¬ 
imately  as  inverse  two  power.  Substituting  an  average  line  of  this  slope  in  Eq.  (37)  reduces  it  to 

A  Kw  =  (1.5  -  /Ox  6.3  (MPasfiH)  (38) 

This  suggests  a  trend  of  A  KTH  from  about  9  to  3  MPa  yfm  as  R  varies  from  zero  to  plus  one, 

independent  of  yield  strength.  This  is  in  rough  accord  with  literature  data  [39,40]. 

The  effect  on  Kk,  however,  is  not  independent  of  yield  strength.  The  model  equates  Kk 
to  a  product  of  the  tensile  instability  strain  and  root  of  process  zone  size.  In  mild  steels,  the 
former  tends  to  vary  inversely  with  yield  strength  level  and,  as  noted  above,  the  latter  also. 
From  this,  one  would  expect  Kk  to  decrease  markedly  as  yield  strength  increases.  This  is  com¬ 
mon  knowledge,  for  example,  the  general  character  of  trend  lines  of  the  Pellini  ratio  analysis 
diagram  [42],  However,  this  argument,  via  association  of  dT  with  grain  size,  also  implies  a  Kk 
increasing  with  grain  size.  This  is  contrary  the  trends  found  by  Stonesifer  and  Armstrong  [42], 
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using  prior  austenite  grain  size  in  a  pressure  vessel  steel.  Perhaps  using  the  grain  size  after  heat 
treatment  would  reconcile  this  inconsistency? 

6.6.  Strain  Rate  Effects 

Inasmuch  as  testing  strain  rate  affects  the  stress-strain  curve,  it  also  influences  growth 
rate  factors  of  the  TLIM.  For  best  record  precision,  our  tests  are  run  at  a  slow  head  speed,  giv¬ 
ing  a  constant-load  strain  rate  of  about  10“ 3  sec-1;  test-fixture  compliance  markedly  retards  this 
in  regions  of  high  tangent-modulus.  With  this,  a  5%  strain  excursion  would  require  about  102 
seconds.  Typical  cycle  durations  in  fatigue  tests  are  a  factor  of  10  to  1000  shorter  than  this. 

Increased  strain  rate  increases  the  elastic  part  of  the  flow  stress.  It  will  also  affect  the 
strain  hardening  rate.  After  Beeuwkes  [43],  body-centered-cubic  steels  will  tend  to  be  9- 
invarient  with  strain  rate.  At  strains  near  the  proportional  limit  changes  in  strain  and  stress  will 
tend  to  cancel,  leaving  the  growth  rate  factor  unaffected.  At  higher  strains,  the  augmented 
stress  will  tend  to  increase  GRF  more  than  the  elastic  strain  augmentation  of  total  strain,  there 
of  diminished  percentage,  will  decrease  it.  However,  the  A-strain  will  be  increased  by  the  elas¬ 
tic  strain  augmentation,  particularly  on  a  log  scale,  in  the  low  AX'  region.  The  net  effect  is  to 
shift  the  da/dN  prediction  to  the  right  at  low  to  intermediate  AK  regions  but  to  merge  all  fre¬ 
quencies  at  higher  AK  regions.  The  air  fatigue  data  of  Vosikovsky  on  X65  steel  [11]  indicates 
such  a  trend  with  cyclic  frequency. 

In  the  hexagonal  titanium  alloys,  Beeuwkes  sees  a  more  moderate  9  dependence  on  strain 
rate  relative  to  <r.  This  would  tend  to  reduce  the  augmentation  of  growth  rate  factor  at  mid-  to 
high  A K  levels,  producing  a  more  continuous  rightward  shift  of  higher  frequency  data.  Recent 
results  of  Chakrabortty  et  al  [44]  indicate  such  an  effect. 
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In  the  face  centered  cubic  aluminum  alloys,  the  strain  rate  sensitivity  tends  to  be  so  low 
as  to  have  little  effect  of  this  sort. 

7.  CONCLUSION 

One  would  like  to  believe  that  the  TLIM,  as  presently  constituted,  is  ready  for  use.  A 
number  of  "pros  and  cons"  weigh  in  evaluation  of  its  utility.  On  the  positive  side:  it  is  firmly 
based  on  measured  stress-strain  data;  it  employs  simple,  limiting-case  assumptions  and 
hypotheses;  it  fits  with  parameters  of  minimum  number,  yet  of  physical  implication;  it 
embraces  all  aspects  of  FCG  and  CFCG  from  threshold  to  instability,  and  all  environments 
including  air  and  vacuum.  On  the  negative  side  it  is  true  that:  the  extra  test  for  stress  strain 
curves  is  burdensome;  the  analysis  is  primitive,  verging  on  mere  dimensional  analysis;  the 
fitting  papameters  are  largely  unverified  by  direct  measurement;  and  in  some  cases,  the  TLIM 
fit  is  imperfect. 

It  is  suggested  that  utility  be  judged  in  two  separate  applications.  As  a  normalizing  pro¬ 
cedure  for  organizing  FCGR  data  with  a  minimum  parameter  system,  it  would  seem  ready.  In 
this  use  it  should  serve  more  effectively  to  focus  the  large  ongoing  effort  in  this  field  by 
providing  a  unifying  basis  for  data  comparison.  The  other  application  is  in  fatigue  life  predic¬ 
tion.  Crack  growth  rate  curves  of  maximum  accuracy  are  needed  for  this,  particularly  in  the 
threshold  region.  Here  it  would  appear  that  in  certain  classes  of  structural  materials,  TLIM  is 
of  useful  accuracy,  hence  could  be  employed.  In  others,  more  research  is  needed  to  determine 
whether  this  is  feasible. 
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